SINS - Small Ionized and Neutral Structures in the Diffuse Interstellar Medium 
ASP Conference Series, Vol. ** VOLUME**, 2006 
M. Haverkorn & W. M. Goss 



Parsec-scale Constraints on the Ionized Interstellar 
Medium with the Terzan 5 Pulsars 



Scott M. Ransom 

National Radio Astronomy Observatory, 520 Edgemont Rd., 
Charlottesville, VA 22901 

Abstract. Over the past two years, we used a series of GBT observations to 
uncover at least 33 millisecond pulsars in the globular cluster Terzan 5 located 
in the Galactic bulge. We now have 32 timing solutions for the pulsars which 
give us precise positions and dispersion measures (DMs) and indicate that the 
DMs are dominated by variations in the integrated electron density along the 
slightly different sight lines towards the pulsars. At a distance of '^8.7 kpc, 
angular separations between the pulsars range from 0.4"— 100" and correspond 
to projected physical separations of 0.01— 4pc, giving us a unique probe into 
the ionized ISM properties on these scales. Our measurement of the DM struc- 
ture function toward Terzan 5 is not inconsistent with Kolmogorov-like electron 
density fluctuations in the ISM on scales ranging from at least 0.2— 2 pc. 



1. Introduction 



Terzan 5 is a rich globular cluster (GC) located in t he Galactic bulge (?=3?8, 
6=1? 7) at a distance (i=8.7±2kpc ( Cohn et"alll2002l ). Deep and repeated ob- 



servations of Terzan 5 using the Green Bank Telescope (GBT) at 2 GHz since 
mid-2004 have uncovered 30 new millisecond pulsars (MSP s) , bringing the total 
known in the clus ter to 33, by far the most of any GC ( Ransom et al.l [20051 : 



Hessels et al] 120061 ). Many individual GC MSPs are in interesting exotic binary 



systems, yet ensembles of MSPs in the same cluster can provide unique science 
as well. Since the pulsars are at identical (within ~0.1%) and known distances, 
and are separated by small angles on t he sky {<iX'—2') , they are unique probes 
of pulsar lumi nosities and masses (e.g. lAndersonlll992 ^). clust er dynamics (e.g . 
Phinnevlll992l ). and ionized gas — both within the cluster dFreire et al.l l200ll ) 



and in the intervening interstellar medium (ISM; lAnderson 19921 ). 



We have recently determined timing solutions for 32 of the 33 pulsars in 
Terzan 5. These solutions provide (among other things), positions on the sky 
(errors are typically ~0.0l" in right ascension and ~0.l"— 0.4" in declination; see 
Figure 1), dispersion measures (DMs; the integrated number of electrons along 
the line-of-sight to the pulsar) with relative errors ^O.Olpccm"^ and absolute 

errors ;^0.1pccm~^, and the apparent spin period derivative. Fobs- The latter 
values are usually dominated by the gravitational acceleration of the cluster 
itself, such that values of Fobs less than or grea ter than zero i ndicate positions 
behind or in front of the cluster respectively ( Phinnev 19921 ). In this work. 



we use the accurate positions and DMs to compute the DM structure function, 
-Cdm('^^)) over a range of angles 1" ^ SO ^ 100", and to constrain the intervening 
ionized ISM on linear scales from 10^ <^ I ^ 10^^ cm. 
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Figure 1. Positions and Dispersion Measures (DMs) of the 32 Terzan 5 
pulsars with timing solutions. Positional errors are typically ~0.01" in right 
ascension, and ~0. 1—0.4" in declination. Relative errors on the DMs are 
^0.01 pc cni~^, while the absolute DMs are known to ^0.1 pc cm^^. The grey 
circle represents the cluster core radius (rc=7.9"). The two pulsars Terzan 5 O 
and Terzan 5 R are separated by only 1.2" on the sky, yet their DMs are 
different by l.lpccm^'^. 



2. The DM Structure Function 

Variations in ISM electron density, Jng , have been sliown to be consistent with a 
Kolmogorov-like wavenumber spectrum of the form PsueiQ) = C^q'^, with (3 '-^ 
11/3 over a very wide r a,nge of wavenumber q = 2Tr/l (corresp onding to scales 
/ from 10*^ — 10^^ cm; e.g. Armstrong. Rickett. fc Spangleilll995l ). Measurements 



of pulsar scintillation, scattering, and DM time variations have been important 
throughout this range in I. 

For simultaneous observations of multiple pulsars separated by small angles 
59 <C Irad (i.e. GC MSPs), and assuming isotropic electron yaria.t ions, the 
electromagnetic phase structure function (e.g. ICordes Sz RickettI [l998l ) depends 
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only on the separation of the sources in the plane of the sky, 6r = d59, and can 
be written ^ 

D^dr) = {rAfUddOr Cl{s) ds, 

where A is the radio frequency, re is the classical electron radius, a = f3 — 2 and 
fa is a constant (when a=5/3, /q,=88.3). Following Cordes (private comm), if 
we assume that the are constant out to distance dgg, and recognize C^d^s 
as the scattering measure SM(dcfr), 

D^{6r) = {vA?^ f-^y SM(4ff)" 
a + 1 Vaeff/ 

For GC MSPs, we observe DM variations as a function of 56 rather than 
phase variations. But we can relate the DM variations to phase variations using 
(5DM = 6(f) /\rf^. Therefore, we effectively measure the phase structure function 
by determining the squared DM differences as a function of angular separation 
59, which is the DM structure function: 



Dj^mm = f Ay SM(deff)^ 

a + 1 Vaeflf/ 



For Terzan 5 (or other sources at large distances in the Galactic plane within 
the ~1 kpc scale height of rig), d ~ (igfj . This means that a good measurement of 
D-om{50) over a range of 56 will determine the power-law index of the electron 
density variations (from its slope a = /3 — 2 on a log-log plot) and an estimate 
of the scattering measure and (from the y-intercept). To be useful for ISM 
constraints, though, measured DM variations must be dominated by the inter- 
vening ISM rather than an ion ized intracluster medium as has been detected in 
47 Tucanae (iFreire et alJ 1200111. 



As noted recently bv iFreire et al.l ( 20051 1. the total DM spread, ADM, of 
the pulsars in GCs with large average DMs (^50 pccm~^) scales roughly linearly 
with DM. If DM variations were dominated by intracluster gas, ADM would be 
roughly independent of distance and the average DM. However, Terzan 5 has 
one of the largest average DMs (~239pccm~^) and the largest ADM (almost 
lOpccm^'^) of any known cluster. If intracluster gas was the caus e, its density 
woul d be a factor 10—20 greater than that found in 47 Tucanae ( Freire et al.l 



20051 ). Another striking indication of the dominance of the ISM contribution to 
ADM over an intracluster gas contribution comes from the pulsars Terzan 5 O 
and Terzan 5 R (see Figure 1). These pulsars are separated by only 1.2" on 
the sky, yet their DMs are different by l.lpccm"^ in the "wrong" way to be 
explained by gas within the cluster. Terzan 5 O, which has the smaller DM, is 
located behind the cluster (as determined by its negative Pobs)^ while Terzan 5 R 
has the larger DM and is likely in front of the cluster. Such a situation can only 
occur if the DM variations are dominated by an intervening ISM that is highly 
structured on linear scales between 0.01— Ipc. 



3. Data Analysis 



With n = 32 independent timing solutions for the MSPs in Terzan 5, we made 
n(n — l)/2 = 496 measurements of ((5DM)^, one for each pair of MSPs (see 
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Figure 2. The DM structure function, Dum{S0), for the 32 Terzan 5 MSPs 
with timing solutions. The 496 measurements of {6DM)^ and their errors are 
the thin crosses. The grey histogram show the number of measurements in 
40 evenly spaced bins. For bins with ^6 measurements, the ((JDM)^) are 
shown as black boxes with error bars. The thick black line shows the linear fit 
to the bins with ^20 measurements, which determined the power-law index 
a = 1.38 ± 0.26. The grey dotted line shows the best fit model assuming a 
Kolmogorov spectrum of electron density variations with a = 5/3. 



Figure 2). To determine Z)dm('5^), we binned the full range of separations, 
projected to the assumed (i=8.7kpc distanc^B to Terzan 5 using d66, into 40 
logarithmically spaced bins. If there were at least 6 measurements of (5DM)^ 
within a bin, we computed the mean-squared DM difference (((5DM)^). Errors 
on (((5DM)^) were estimated using the measured DM distribution for Terzan 5 
and the number of measurements per bin, according to normal sampling theory. 

For bins with >20 measurements of (JDM)^, we fit a simple linear (in 
log-log space) model to the data: log^o ((((^DM)^)) = alogiQ{d69) + b. For 
the first model, we allowed both a and b to vary and measured a=1.38ib0.26. 
Since the measured power-law index of the turbulence is only ~l-cr away from 
the value for Kolmogorov fluctuations, we determined b by assuming q=5/3, 
and re-fitting to get 6=— 29.893±0.030. Applying appropriate unit conversions. 



^ Since we are probing the intervening ISM, an argument can be made to use a fraction of the 
distance to Terzan 5 rather than the total distance, d. 
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SM = lo(''+29-i66) ^ 0.188±0.013kpcm-20/3, and therefore = SM/d^s 
0.022 m-20/3. 



Discussion and Conclusions 

It is interesting to compare our estimate of tlie SM to wards Terzan 5, . 188 ± 
0.013 kpcm"^'^/'^, to values determined in other ways. Nice &: Thorsett ( 19921 ) 



measured a pulse broadening timescale due to scattering, Tp, for pulsar Terzan 5 A 
at 685 MHz of 0.70 ms. This can be converted to a pulse broadening based SM 
estimate, SM,-, using r„ = 0.90 ms SMr^^z^~^^''^d, with u in GHz and d in kpc 



(|Cordesl[2003 ). which gives SMT-=0.033kpcm "^^Z^, This value is smaller by a 



factor of ~6 than our measurement, although we should expect some systematic 
variation due to the different weighting factors used by the methods over what is 
almost certain ly an inhomogeneous electron distribution along the line-of-sight 
( Corded [2OO2I ). It could also be that there is an upturn in the power spectrum 



between the small scales probed by Tp and the much longer scales relevant to 
the DM structure function. 

Similarly, we can use either the known (DM) ~ 239pccm~3 or distance to 
Terzan 5, (i=8.7kpc, as inputs to the NE2001 Galactic electron density model, 
which can estimate SM (|Cordes fc Lazioll2002l l . If we specify d, NE2001 predicts 



SMrf=0.79kpcm-20/3, while if we specify DM, SMDM=0.29kpcm-2°/3^ ^he 
latter value is quite close to our measurement, and significantly closer than the 
estimate based on the distance to the clusteo 

It is quite interesting that the measured power-law index of the DM struc- 
ture function for linear scales ranging from ~0.2— 2 pc (or wavenumbers, 10~^^ < 
q < 10^^^ m^^) is not inconsistent with Kolmogorov fluctuations of the ISM. 
These scales are 10^ — 10^ times longer than those probed by most other pul- 
sar techniques. When we use our estimates of and assume /3 = 11/3 
to determine Psn^{q), the resulting values are amazingly close to extrapola- 
tions based on scintilla tion and other measur ements at wavenumbers between 
10~^3 < q < 10~^ m~^ ( Armstrong et al. 19951 ). implying that Kolmogorov-like 



density fluctuations are present in the ISM over size scales covering almost ten 
decades. It is intriguing to note the possible downturn of Ddm{S6) at linear 
scales of 2— 3pc. While this is likely due to the poor statistics at the large-angle 
end of the structure function, it could also indicate saturation of the structure 
function upon reaching the outer scale (i.e. the largest size scale where power is 
input into the turbulent process governing the ISM). 

Given the large number of pairwise measurements provided by the Terzan 5 
MSPs, one of our next tasks will be to check for evidence of anisotropic electron 
density fluctuations over the same size scales. In addition, ongoing searches of 
GCs using the GET and other telescopes have resulted in several other bulge 
GCs with substantial numbers of pulsars where these techniques might be useful 
in the future (for instance, M28 currently has 11 known pulsars). 



^If the distance to Terzan 5 is really 8.7 kpc, NE2001 greatly overestimates the integrated 
electron density towards the cluster. For that distance, NE2001 predicts DM~530 pc cm~'^. 
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